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I. INTRODUCTION 



Finite impulse response (FIR) filters have seen widespread use in the field of 
two-dimensional digital signal processing. Reasons for this include the inherent 
stability of FIR filters, the ability to achieve a linear (or zero) phase characteristic 
in the frequency response, and relative ease of design. Traditional approaches to 
the design of such filters include the window method, in which a smoothing window 
is applied to the Fourier series coefficients of the ideal frequency response; the fre- 
quency transformation method, in which a one-dimensional (1-D) prototype filter 
response is mapped into a function of two frequency variables using an appropri- 
ate transformation function; and the frequency-sampling method, in which desired 
frequency response values are specified at certain sample points in the frequency 
domain. Of the three methods described above, the frequency-sampling method is 
probably the least utilized and the least understood, and hence warrants further 
investigation. This thesis deals with a comparative analysis of various methods of 
designing two-dimensional (2-D) FIR filters using frequency-sampling techniques. 
In the analysis, some new methods for frequency-sampling design are presented. 

The 2-D frequency-sampling design method is, in essence, an interpolation 
problem in two variables where the two variables are the two frequencies of interest, 
denoted as uj x and u> 2 . The frequency response is specified at selected sample 
points in the (w!,w 2 ) plane, and is expressed in terms of linearly independent 
basis functions. The solution to this interpolation problem involves finding the 
coefficients associated with the respective basis functions in order to meet the 
interpolating conditions. These coefficients are, in general, closely related to the 
impulse response. The various approaches to the frequency-sampling design 
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problem involve differing constraints placed on sample location in the (uq , w 2 ) plane 
and differing sets of basis functions. Because the resulting frequency response takes 
on exactly the prescribed value at each frequency sample location, the response is 
somewhat predictable at the onset of the design process. A disadvantage of 2- 
D frequency-sampling design, however, is that the existence and uniqueness of a 
solution is dependent upon the location of selected frequency samples and the basis 
functions specified; in other words, degenerate cases may arise in which no solution 
exists to the bivariate interpolation problem. 

FIR filter design techniques denoted as “frequency-sampling” traditionally in- 
volve sampling a desired frequency response at equally-spaced intervals and then 
solving for the impulse response coefficients using an inverse discrete Fourier trans- 
form (IDFT) [Ref. l] . For the 2-D case, the process is similar, with the desired 
frequency response sampled at equally spaced points on a Cartesian grid in the 
(uq,u> 2 ) plane [Refs. 2, 3]. This approach will be referred to as the “uniform 
sampling” method in this thesis. There are two significant benefits to such an ap- 
proach in two dimensions. First, with frequency samples placed at the vertices of 
an N x N Cartesian grid, the existence of a unique filter satisfying the interpolating 
conditions and possessing an impulse response with an N x N region of support is 
guaranteed. Second, use of the IDFT, or alternatively, the fast Fourier transform 
(FFT) algorithm, assures a computationally efficient design process. The princi- 
pal disadvantage of the approach is that there is no flexibility in the placement of 
frequency samples. The approach is presented in more detail in Chapter II. 

A different approach to frequency-sampling design of 2-D FIR filters involves 
allowing frequency samples to occur anywhere in the (uq,w 2 ) plane and fixing 
the region of support for the filter impulse response. This approach is presented 
in Chapter III and will be denoted as the “arbitrary sampling” technique. The 
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method is quite flexible since samples can be arranged to enhance certain desired 
features in the resulting frequency response, such as flat passband, sharp roll-off, 
equi-ripple behavior, etc. The design, however, involves the solution of a large 
system of linear equations for a filter of any practical order and hence it is quite 
computationally intensive. Additionally, with no constraints on the location of 
samples, degenerate cases may arise in which no solution exists which meets all 
interpolating conditions for the specified impulse response region of support. 

A novel 2-D frequency-sampling design technique is presented in Chapter IV, 
and is referred to as the “nonuniform sampling” method. This method places 
certain constraints on frequency sample location but still possesses much greater 
flexibility than the uniform sampling approach. Rather than solving one large 
system of equations, this technique involves the solution of several smaller systems 
of linear equations. Given the constraints on the sampling geometry, existence and 
uniqueness of a solution is guaranteed. Additionally, the design method is less 
computationally intensive than the arbitrary sampling case but more intensive, in 
general, than the uniform sampling case. 

In comparing and contrasting these various approaches to 2-D frequency- 
sampling design, several interesting points arise. First, frequency-sampling design 
techniques involve some sort of decision-making process in determining where to 
place the samples in order to best match the desired frequency response. Methods 
which are more constrained, such as the uniform sampling approach, involve a lesser 
degree of decision-making and hence can be considered easier to use from that 
perspective. Second, by reducing the number of constraints on the placement of 
frequency samples, flexibility in matching the desired frequency response is gained. 
Finally, and perhaps most significantly, computational efficiency in determining 
filter parameters can be directly related to the number of constraints on sample 
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location. Of the techniques discussed here, those which are more constrained tend 
to be more computationally efficient. Thus there are tradeoffs to be considered 
when selecting one of these techniques. In this thesis, a detailed analysis of these 
concepts is performed and the quality of these frequency-sampling methods in 
matching desired filter characteristics is investigated through design examples. 
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II. DESIGNS WITH UNIFORM SAMPLING ON A CARTESIAN GRID 



The approach to frequency-sampling design of 2-D FIR filters discussed in this 
chapter is the one referred to as “frequency-sampling” in various texts [Refs. 1, 

2]. It will be referred to the “uniform sampling” method through the remainder of 
this thesis. This method has significant computational advantages through use of 
an inverse discrete Fourier transform (IDFT) approach. Additionally, this method 
will never result in any degenerate cases, i.e. situations in which a filter satisfying 
each of the interpolating conditions cannot be found. However, due to constraints 
on the location of samples, the method is rather inflexible. 

A. APPROACH 

The transfer function of a two-dimensional FIR filter with impulse response 
support in the first quadrant can be defined as 

N i-l N j — 1 

H(z l ,z 2 ) = ^2 h(n 1 ,n 2 )z~ ni z~ n \ (2.1) 

n i = 0 n 2 = 0 

The frequency response of a such a filter is obtained by evaluating the transfer 
function on the unit bicircle through the substitutions z x = e ,Wl and z 2 = e ;W3 , 
yielding 

N i-l N,-l 

H(u 1 ,u 2 )= ^2 h(n 1 ,n 2 )e _ywini e _iW5n5 . (2.2) 

fii=0 n 3 = 0 

The frequency-sampling filter design problem involves determination of a sequence 
h{n i , n 2 ) which will result in a filter matching a desired frequency response H d (u^ , w 2 ) 
at specified sample points in the (uq,a; 2 ) plane. 

The initial step in the design process is to choose the filter order, and hence 
the parameters N x and N 2 . Often it is desirable to specify Ni = N 2 = N to 
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ensure a result with equal dependence on the frequencies uq and u > 2 . The desired 
frequency response i/ (i (uq ,u> 2 ) is sampled at the vertices of a uniformly spaced 
Ni x N 2 Cartesian grid on the region {0 < uq < 27r; 0 < ui 2 < 2x} of the (uq,tu 2 ) 
plane. The resulting frequency samples can be expressed in 2-D sequence form as 



H (ki , k 2 ) — H d (uq , ) | 3«*, , 

l ~ N i • N i 

Aq =0,1 ,Ni — 1 and k 2 = 0, 1, . . . , N 2 — 1. (2.3) 

The resulting sampling geometry is illustrated in Figure 2.1. 

Wo 




0 2 rr 4 n 2 n ( N 1 - 1 ) 

* i N , N , 

Figure 2.1. Arrangement of Samples for Uniform Sampling Method. 



Because an IDFT approach is to be used in computing /i(n 1 ,n 2 ), and since the 
IDFT is defined for first-quadrant sequences only, the desired frequency response 
should possess a linear phase characteristic such that 

(2.4) 
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where 



My = 



Ny/2 N v even 
(Ny — l)/2 Ny odd 



( 2 . 5 ) 



M _ f N 2 /2 N 2 even 

2 \ (N 2 — l)/2 N 2 odd ' 

Additionally, H d (u>y, u> 2 ) must possess the symmetry conditions 

\H{ky,k 2 ) = \H{ky,N 2 -k 2 )\ 

= | H(Ny~ky,k 2 )\ 



( 2 . 6 ) 



( 2 . 7 ) 



in order to ensure a purely real impulse response. 

The next step is to compute h(ny,n 2 ) by taking the IDFT of the sampled 
desired frequency response: 



h{ny,n 2 ) = IDFT [H{ky,k 2 )\ 



1 

nJT 2 



W, - 1 ATj-l 



k x =0 fc 3 = 0 



( 2 . 8 ) 



where 

W Nl =e~ 3 "T and W N , = e _J *T. 



The resulting sequence h(riy ,n 2 ) is purely real and represents the impulse response 
of the designed filter. Such a filter possesses a linear phase characteristic. A zero- 
phase characteristic can be obtained by shifting the impulse response so that it is 
centered about the origin such that 



ti(my,m 2 ) = h(my + My,m 2 + M 2 ), 



( 2 . 9 ) 



where h 1 (m 1 ,m 2 ) represents the zero-phase impulse response. 
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To see the relationship between the resulting filter frequency response and the 
specified sample sequence H(k 1 ,k 2 ), Eq. (2.8) is substituted into Eq. (2.2): 

Wl-lfT)-lr tfi-Ufi-1 

*(«..<*)= E E wir'E E 



NiN 2 

m = 0 n 7 = o L k 1= = 0 k 7 = 0 



e ~ joj i n l e -ju 7 n 7 



N 7 - 1 



= 7fV E E E; E (W^-e-'-r 



^iV 2 



fc i = 0 fc 2 =0 



fii = 0 



n 2 = 0 



iNTi — 1 N 7 - 1 

= E E 

Ac. =0 *J = 0 



_ p-iuiNi 



) 



H(k t ,k 2 ) 



N x N 2 






The above expression can be rewritten in the form 

tf.-iWj-i _ 2 

H(u ,,wj) = E E *(*.. **)*(■». 

fci =0 fc 2 =0 1 

where 

(1 - e“ iWlWl )(l - e - - 7 "’ ) 



(2.10) 



( 2 . 11 ) 






./Vx iV 2 (l — e _ -’ w ‘)(l — e - - 7 "’) 



which simplifies to 



$(u>i,u> 2 ) 



sin^jiVj /2) sin(w 2 jV 2 /2) 



iVj N 2 sin(u; 1 /2)sin(u; 2 /2) (2.12) 

. g -iwi( iLi 5 rJ ’) e -- 7U ’5( iL ¥ :i )e J ’ r,Cl(1_ «r ) e i7rfc,(1- "T 1 

for the case where N x and N 2 are both odd integers. This analysis is a 2-D 
extension of operations performed for the one-dimensional case [Ref. 1: pp. 97- 
98]. Since 

&(*lk —k\ = { 1 (*i,A: a ) = ( 0 , 0 ) 

K N, " N 2 2> \0 fc, = 1,2,...,^ - 1 and k 2 = 1, 2, . . . , N 2 - 1 ’ 

the frequency response is matched exactly at the original sample points. That is, 



H (u>x , w 2 ) | _ 1 _ a<r * 3 — H (fci ,k 2 ), 

»1 ’ 3_ "3 



(2.13) 
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k 1 = 0, 1, . . . , N 1 - 1 and k 7 = 0, 1, . . . , N 2 — 1. 



Thus, the design process entails a bivariate interpolation in the variables u> l and 
u> 2 , and, due to use of the IDFT algorithm, the interpolating functions are of the 
form of Eq. (2.12). 

In using this filter design method, the choice of the desired frequency response 
H d (u x ,cj 2 ) has significant effects on the quality of the resulting filter. For instance, 
if an “ideal” lowpass filter characteristic with the circular symmetry shown in 
Figure 2.2 is used for H d (u 1 ,u> 2 ), the passband and stopband of the resulting filter 
frequency response exhibit substantial ripple. 




Figure 2.2. Ideal Lowpass Filter Frequency Response. 
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Figure 2.3 illustrates contour and perspective plots of the frequency response 
arising from such a design. In this case, a 15 x 15 point filter was designed by 
uniformly sampling the desired frequency response of Figure 2.2. 

The frequency response of the resulting filter can be improved considerably if 
a transition region is specified for H d {u> l ,u> 2 ). Increasing the width of such a tran- 
sition region reduces passband and stopband ripple but does so at the expense of a 
less sharp roll-off characteristic in the resulting filter response. The values assigned 
to H(ki,k 2 ) for points (k x ,k 2 ) located in this transition region significantly influ- 
ence the behavior of the resulting filter. The selection of these values in order to 
satisfy some specific optimization criterion for the frequency response, such as best 
equi-ripple approximation to an ideal response, is not trivial. Since H(u> i,u/ 2 ) is a 
linear function of the values H(k 1 ,k 2 ), as seen in Eq. (2.10), linear optimization 
techniques can be employed in selecting the transition region values of H(k 1 ,k 2 ). 
Such techniques for the one-dimensional case are rather involved [Ref. 4]; for the 
2-D case, the problem is much more computationally intensive [Ref. 2], and hence, 
has not been addressed. Specifying H d (uq , u 2 ) to have a smooth, linear transition 
between the passband and stopband regions is the simplest thing to do and gen- 
erally works well. The magnitude characteristic of a desired frequency response 
with a linear, circularly symmetric transition region is shown in Figure 2.4. A fil- 
ter can then be designed by sampling this response on the uniform Cartesian grid 
specified in Eq. (2.3). Contour and perspective plots of a 15 x 15 filter designed 
in this manner are shown in Figure 2.5. In comparing this result to that of Figure 
2.3, it is evident that inclusion of a transition region greatly reduces passband and 
stopband ripple. 
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(a) 




(b) 



Figure 2.3. Frequency Response Arising From Sampling Ideal Lowpass 
Characteristic With No Transition Band. 
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H (cjI ,co2) 



1.0 




Figure 2.4. Desired Lowpass Frequency Response Characteristic 
With Linear Transition Band. 



B. EXISTENCE AND UNIQUENESS 

A significant property of the uniform sampling design is that a unique filter 
satisfying each of the interpolating conditions specified in Eq. (2.3) is guaranteed 
to exist. In general, existence and uniqueness of such an interpolation result in two 
dimensions is not guaranteed; this can be a significant theoretical and practical 
problem. The point is discussed in further detail in Chapter III. A proof of existence 
of a unique h(n 1 ,n 2 ) using the frequency-sampling method with uniformly spaced 
samples follows. 
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(a) 




(b) 



Figure 2.5. Frequency Response Arising From Sampling Desired Lowpass 
Characteristic With Linear Slope Transition Band. 
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Equation (2.8) is used to determine h(n 1 ,n 2 ). This can be rewritten as 



^,-lr W,-l 



Y, ^ 



N- 



k,= 0 L fc, = 0 



rr^-fcini 
rr W, 



(2.14) 



Now define 



1 Ni-l 

G(k„n,) H(k,,k 2 )W^”". (2.15) 

2 fc , =0 

Substituting Eq. (2.15) into Eq. (2.14), the following expression for the impulse 
response results: 

N i ~ 1 

h(ni,n 2 ) = — G(k 1 ,n 2 )W~* ini . (2.16) 

1 k, = 0 

Thus the 2-D IDFT of H (Aq , k 2 ) can be determined in a two-step process involving 
one-dimensional transforms known as row-column decomposition [Ref. 5: pp. 75- 
76]. 



The first step in the row-column decomposition process is to perform the 
solution of a total of N x systems of N 2 equations of the form of Eq. (2.15) using a 
1-D IDFT algorithm. In matrix form, this can be written as 



G(Ah, 0 ) I 




r 

N-, 


w°„ 

N i 


••• K, 1 




r Hik^o) i 


G(Ah,l) 


1 


w° 

N-, 


wal 


... IT" 1 "’" 1 ' 

ATq 




H{k lt 1 ) 


-G{k 1 ,N 2 - 1 ) - 




W 9 . 

N, 




w - 1 




.H{k lt N 2 - 1 ). 

( 2.17 



which can be expressed in vector-matrix notation as 



G = ~W 2 H. (2.18) 

It can be shown [Ref. 6: pp. 44-45] that a scaled version of W 2 , namely W 2 / \J iV 2 , 
is unitary and hence non-singular. Therefore a unique solution for G exists for each 
k<, i = 0, 1, . . . ,Ni - 1. 
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The second step in the computation involves the solution of N 2 systems of N 1 



equations of the form of Eq. (2.16). In matrix form, this becomes 



h(0,n 2 ) 




rw° 

N , 


W° N 

N ! 


W° 1 

Ni 




- G(0,n 2 ) * 


h(l,n 2 ) 


1 


W° 

yy N t 


w~ l 

yy Nl 


ry-(Wi-i) 

yV N l 




G(l,n 2 ) 


.h(Ni — 1 , n 2 ). 


“ N~i 


w° 

lvy n i 


Vy N l 


N i J 




. G(Ni — 1, n 2 ) . 

(2.1 1 



which in vector-matrix form is 



h = -^-WjG. (2.20) 

Ni 

Since G was uniquely determined for each value of k x , and since W t is also 
non-singular, a unique solution for h exists for each ny, j = 0, 1, . . . , iV 2 — 1. 
Therefore, a unique solution for h(n 1 ,n 2 ) is obtained from the sequence H(k 1 ,k 2 ), 
which represents the samples of a desired frequency response at the vertices of a 
uniform Cartesian grid in the (wi,w 2 ) plane. 

C. COMPUTATIONS 

While there is no flexibility in the choice of sample point locations using this 
IDFT approach, the method does offer some significant computational advantages. 
One criterion often used to measure computational efficiency in solving a problem 
is to determine the number of complex operations involved. While there exist many 
methods of computing the 2-D IDFT, the use of row-column decomposition and 
the IDFT directly (vice the FFT algorithm) is a reasonably efficient approach for 
the design of filters of practical order and is used in the following analysis. 

The initial step in the row-column decomposition approach is to solve a total 
of Ni systems of the form of Eq. (2.17). Solution of each system via the IDFT 
involves a total of iV| multiplications and additions. Therefore, a total of A^iVl 
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floating point operations are required in this step, where a floating point operation 
is considered one multiplication and one addition. 

The second part of the row-column decomposition approach involves the so- 
lution of N 2 systems of the form of Eq. (2.19). Solution of each system via the 
IDFT requires iVf operations. Therefore, N 2 floating operations are involved 
in the second step. 

The total number of floating point operations required in the solution for the 
impulse response coefficients of an Ni X N 2 frequency-sampling FIR filter design 
using uniform sampling on a Cartesian grid and a row-column decomposition, IDFT 
approach, is 

C uniform = N 1 N 2 +N 2 1 N 2 . (2.21) 

For the case in which N x = N 2 = N, this reduces to 

C UIllform =2 N\ (2.22) 

In summary, the uniform sampling, IDFT approach to 2-D frequency sampling 
filter design has the benefit of computational efficiency, yet lacks flexibility in the 
choice of frequency samples. The next two chapters will present methods which 
achieve greater flexibility by allowing nonuniform spacing of samples. 
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III. DESIGNS WITH ARBITRARILY PLACED SAMPLES 



In Chapter II a frequency-sampling 2-D FIR filter design method was pre- 
sented in which samples were placed on a uniform Cartesian grid in the (uq , u 2 ) 
plane. Due to the constraints on sample placement, however, the method was not 
particularly flexible. If the desired frequency response H d (uq , uq) was fixed, no free 
parameters were available in optimizing the frequency response of the resulting fil- 
ter. In this chapter, an alternate frequency-sampling design technique for 2-D FIR 
filters is presented in which samples are no longer constrained to a uniform Carte- 
sian grid but can be placed anywhere in the (uq , u 2 ) plane. This approach will 
be referred to as the “arbitrary-sampling” method. Allowing the samples to occur 
anywhere in the (uq , uq) plane introduces additional free parameters; hence, linear 
optimization techniques can be employed without varying H d (uq , uq ) , the desired 
frequency response to be sampled. Since the locations of the frequency samples are 
much less constrained, the method provides greater potential for achieving a desir- 
able frequency response characteristic in the resulting filter. Since the placement 
of samples on a uniform Cartesian grid is a special case of the arbitrary-sampling 
method, the arbitrary-sampling method can achieve results at least as good as 
the uniform-sampling method for a fixed impulse response region of support and 
a fixed number of samples. However, the arbitrary-sampling method possesses a 
few significant drawbacks, such as computational complexity and the possibility of 
degenerate cases in which the interpolating conditions cannot be met. These must 
be considered in the selection of a method for filter design. 
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A. APPROACH 



The arbitrary-sampling method can be used to produce 2-D FIR filters with 
a variety of impulse response regions of support. This region of support is one of 
the parameters to be specified at the outset of the design process. For illustration 
purposes, a filter with an N x x N 2 rectangularly-shaped impulse response region 
of support centered on the origin of the (n. l5 n 2 ) plane is specified, where N l and 
N 2 are odd integers. The transfer function of such a filter can be expressed as 

M i Mi 

H{z 1 ,z 2 )= 22 22 h i n i,n 2 )z~ ni z~ n \ (3.1) 

n i = — A T i n 3 = — AT 3 

where h(rii,n 2 ) is the impulse response to be determined, M x — ( N x — l)/2 and 
M 2 = ( N 2 — l)/2. The parameters M x and M 2 will be used often in the remainder 
of this analysis, so the above relationships are worth noting. Evaluation of Eq. 
(3.1) on the unit bicircle yields the following expression for the frequency response: 

M ! M 3 

H{u 1 ,u 2 )= Y, Y, M n i ,n 2 )e~ JUini e~ :u * n \ (3.2) 

n x = — AT 1 n 3 = — AT 3 

To ensure a zero-phase characteristic in the resulting frequency response, 
h(ny,n 2 ) is specified to possess quadrant symmetry such that 

h(u-! , n 2 ) = h(—n 1 , n 2 ) = h(n 1 , —n 2 ). (3.3) 

Application of Eq. (3.3) to Eq. (3.2) produces the following expression: 

Mi 

H(u>y , u> 2 ) — h(0, 0) + '22 2h.(n 1 ,0)cosw 1 n 1 -|- 

n 1 = 1 

AT 3 AT 1 AT 3 

22 2/i(0, n 2 ) cos u 2 n 2 + 22 X/ > n 2 ) cost^ri! coscj 2 ™ 2 - 

n 3 = 1 fli = l «3=1 

(3.4) 
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This can be rewritten as 



M ! 

cj 2 ) - 



M 2 

^ a(n L ,n 2 ) cos uqn! cos u> 2 n 2 , 



where 

a[n i ,n 2 ) 

The basis functions {& (w! , 



ni = 0 n 2 = 0 

{ h(ni,n 2 ) rii = 0,n 2 = 0 
2h(n 1 , n 2 ) rii ± 0, n 2 = 0 
2 h(n 1 ,n 2 ) n x = 0,n 2 ^ 0 
4h(rii ,n 2 ) n x ^ 0, n 2 ^ 0 

w 2 )} are defined such that 



rpi (cj t , u 2 ) = cos <jj x n x cos u> 2 n 2 , 



(3.5) 



(3.6) 



(3.7) 



where 

i = n x [M 2 + 1) + n 2 + 1; n 1 = 0, 1, . . . , M x ; 

(3.8) 

n 2 = 0 , 1 ,..., M 2 . 

By arranging the 2-D sequence a{rii,n 2 ) into the one-dimensional sequence a(i ) , 
the frequency response can then be expressed as 

R 

H(uji,uj 2 ) = ( 3 - 9 ) 

t= 1 

where 

R = [Mi +1)(M 2 +1). (3.10) 

From Eq. (3.9), the elements of the sequence a(i) represent a total of R free 
parameters to be determined. Therefore, the frequency response is to be specified 
at a total of R distinct sample points in the (u; 1 ,cj 2 ) plane. 

The initial step in the design process is to specify H(u> 1 ,u> 2 ) at R distinct 
sample points in the region K = {(t^i,w 2 ) : 0 < u> x < tt; 0 < u 2 < 7r} . An example 
of such a sampling geometry is illustrated in Figure 3.1. 

No constraints other than the distinctness requirement are placed on the loca- 
tion of samples within this region. The form of Eq. (3.5) implies that the frequency 
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Figure 3.1. Possible Arrangement of Samples for Arbitrary Sampling Design. 



response possesses quadrantal symmetry and periodicity of in both frequencies. 
Hence, H(u> 1 ,u> 2 ) is defined for all (uq,w 2 ). 

The next step in the design process is to solve for the sequence a(n! , n 2 ). Appli- 
cation of Eq. (3.9) at each specified sample location {(uq ,u> 2 ), ; * = 1,2,..., 72} 

results in the following system of linear equations: 



■ H(Ui,U 2 )i ' 
,w 2 ) 2 




01 (u»i , W 2 ) 1 02 (^1 5 ^2 )l 

0l(Wl,W 2 ) 2 ^2(Wi,W 2 ) 2 


••• 0fl(Wl,W 2 )i ' 

... 1p R (bJ i,W 2 ) 2 




r«(i)- 

a(2) 


1 

e 

e 

10 

i 




-0l(w lt W 2 )ii 02 (Wi,W 2 )/l 


... 0/j (u, ,U 2 ) R . 




.a(R). 



(3.11) 
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In vector-matrix notation, Eq. (3.11) becomes 



H = Ba, (3.12) 

in which B is of dimension R x R. The sequence a(n 1 ,n 2 ) results from solving the 
relation 

a = B 1 H. (3.13) 

The final step in the process is to compute the filter impulse response h(n 1 , n 2 ) 
using the relationship of Eq. (3.6). The resulting impulse response has an N x x N 2 
region of support. Equation (3.4) suggests a transfer function of the form 

M, Mi 

H{z x ,z 2 ) =h{0,0) + M«1,0 )«* +z; n ')+ J2 Mo ,n 2 ){z^ +z; n ') 

n i = 1 n. 2 = 1 

M, 

+ Hni t rh)(zV +2r" , )«’ + z 2 n '), 

fi!=l n 2 = 1 

(3.14) 

which can be used in implementing a 2-D filter structure with a minimal number 
of gains. 

This same approach can also be employed if other than rectangular impulse 
response regions of support are specified. Again, a system of linear equations of 
the form of Eq. (3.5) is solved, but the limits on the summations are altered to 
reflect the desired region of support for h(n 1 ,n 2 ). 

To summarize the arbitrary-sampling design process, the following steps are 
performed: 

(1) Fix the impulse response region of support, and hence, the filter order. 

(2) Specify the desired frequency response at selected sample points on the 

region K = {(u^Wj) : 0 < < tt; 0 < u> 2 < rr}. The number of specified 

samples must equal the number of independent impulse response values 
(i.e., those which are not fixed by the symmetry conditions). The 
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locations of samples within the region is unconstrained except for 
the condition that each sample location be distinct. 

(3) Solve a system of linear equations to obtain the resulting impulse 
response h(n 1 ,n 2 ). 

B. EXISTENCE AND UNIQUENESS 

There is a fundamental deficiency in the arbitrary-sampling approach to the 2- 
D FIR frequency-sampling filter design problem. For specified frequency response 
samples placed arbitrarily in region K , there is no guarantee that a unique filter 
possessing a number of independent impulse response samples equal to the number 
of interpolating conditions, and which satisfies each of those interpolating condi- 
tions, exists. This is a problem inherent to interpolation in two dimensions which 
does not arise in the one-dimensional case. 

For a unique solution to the system of equations indicated by Eq. (3.12) to 
exist, the basis functions {V’< (^i , cj 2 )} must satisfy the Haar condition. A set of 
vectors in n-space satisfies the Haar condition if every set of n of them is linearly 
independent [Ref. 7: p. 45). Therefore, to satisfy the Haar condition, any set of 
R characteristic vectors defined on R distinct points in region K of the (wj , w 2 ) 
plane must be linearly independent. For the arbitrary-sampling method developed 
in this chapter, a characteristic vector {^((cji , w 2 ),)} associated with a particular 
sample location (c^q , w 2 ),- can be defined as 

®((wi,w 9 )i) = [^((wi.wjjj) V' 2 ((wi,w 2 ) i ) ... if) R ((ui,u 2 )i) ] T . (3.15) 

It is known that for a region K' of the two-dimensional frequency plane, there are no 
nontrivial sets of two-dimensional basis functions which satisfy the Haar condition 
[Ref. 8: p. 494). Therefore, there is no guarantee that the sets of characteristic 
vectors formed by evaluating Eq. (3.15) at arbitrarily selected sample points in 
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region K of the (uq,(j 2 ) plane are linearly independent. Since the square matrix 
B in Eq. (3.12) can be expressed in terms of the characteristic vectors as 

B = [ $(^, 0 ) 2 ) 1 ) ’® r ((wi,w 2 ) 2 ) ••• ®((wi,W2)ji) ] T , (3.16) 

and since the characteristic vectors may be linearly dependent, B may be singular 
and a solution to Eq. (3.13) may not exist. Such a case is termed a degeneracy, 
and illustrates a basic problem associated with the arbitrary-sampling method. As 
can be inferred from the form of the basis functions, the occurrence of a degeneracy 
is dependent solely upon the location of samples; the specified frequency response 
values at the sample points do not affect this. 

A simple example illustrates how a degeneracy can arise. The frequency re- 
sponse of a filter is specified at four sample points in the (cj 1 ,cj 2 ) plane: 

1. iJ(0, 0.67t) = 1 2. H(0An,0) = 1 

3. H( 0.47r,7r) =0 4. /f (7r, 0.67t) = 0 

These sample locations are shown in Figure 3.2. 

The parameters and M 2 are both set equal to 1 so that the number of 
samples in the sequence a(n 1 ,n 2 ) equals the number of frequency samples. Sub- 
stitution into Eq. (3.11) yields 



1- 




"1 COS0.67T 1 COS0.67T 




'g(0, 0) ' 


1 




1 1 COS0.47T COS0.47T 




g(0, 1) 


0 




1 —1 COS 0.47T — COS0.47T 




g(1,0) 


0 . 




. 1 COS 0.67T —1 — COs0.67T_ 




Lg(1,1)J 



Some calculations reveal that the square matrix above has a determinant of zero; 
therefore, it is not invertible and no unique solution for a(n.!,n 2 ) exists. 

To determine whether degenerate cases present a serious shortcoming of this 
method, it is of interest to determine how often degeneracies are likely to arise. To 
investigate this, an APL computer program was developed to count the number of 
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Figure 3.2. Arrangement of Samples for Design Example. 

degenerate cases arising from the arbitrary selection of sample points {(u>x ,cu 2 ), } 
on region K of the (u>!,w 2 ) plane. A total of R samples were arbitrarily placed 
at distinctly different vertices of a grid of specified density on region K for each 
case. The placement of samples was performed in a pseudo-random fashion using 
the APL “roll” function. The density of the grid, in effect, is determined by the 
number of decimal places carried in the frequency values. In the analysis, two 
parameters were varied: the dimensions of the specified filter impulse response, 
and the grid density. For each situation, one hundred separate sets of arbitrary 
samples were investigated. The program then determined the number of cases in 
which the sampling geometry resulted in a degeneracy. Due to the large number of 
computations involved, analysis was limited to filters with square impulse response 
regions of support only. The results are listed in Table 3-1. 
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These results support two concepts. First, for a fixed grid density, as the 
order of the specified filter increases (and hence the number of specified frequency 
samples), the likelihood of the occurrence of a degeneracy increases. Secondly, for a 
fixed filter order, as the density of the grid on which samples are placed increases, 
the likelihood of a degeneracy occurring decreases. For filters of order 17 X 17 
or less, the probability of a degenerate case is minimal provided frequencies are 
specified to at least the nearest 0.05 radians. 

TABLE 3-1 

NO. OF DEGENERACIES OCCURRING IN 100 TEST CASES 





GRID DENSITY (RADIANS) 


FILTER ORDER 


R 


0.2 


0.1 


0.05 


5x5 


9 


5 


1 


0 


9x9 


25 


10 


0 


0 


13 x 13 


49 


17 


1 


0 


17 x 17 


81 


56 


2 


0 



It is significant to note that for the above tests, frequency samples were cho- 
sen in an arbitrary manner. However, when designing a filter to approximate 
a certain desirable frequency response characteristic, the selection of points is 
likely to be done in an unconstrained, but not arbitrary, manner. Certain pat- 
terns in sample selection, such as selecting octant symmetric sample points where 
(w! ,u> 2 ); = (u> 2 >u>i )y, have a greater tendency to introduce linear dependencies 
within the set of characteristic vectors than when sample points are selected at 
random. This concept is difficult to show in a controlled experiment; however, in 
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the course of this research, degeneracies arose when attempting filter designs much 
more often than the arbitrarily placed sample test results would indicate. 

Such a degeneracy did arise in the following design example. In this case, a 
lowpass filter is desired with a circularly symmetric frequency response such that 
it possesses unity gain in the passband region {0 < \J < 0.2n} and zero 
gain in the stopband region {0.47T < < 7r}. The region of support for 

the impulse response is specified to be 17 x 17, implying Af x = M 2 = 8. Therefore, 
a total oi R = (Mi + 1 )(M 2 + 1) = 81 frequency samples are to be specified. Figure 
3.3 shows the location of 81 points at which the desired frequency response is to 
be sampled initially. 




Figure 3.3. Arrangement of Samples for Lowpass Filter 
Design, Case 1. 



32 



Note that the arbitrary-sampling method does not require any samples to be 
placed in the transition (don’t care) region. This allows the concentration of more 
samples in the passband and stopband regions. This particular choice of sample 
locations, however, was found to result in a degeneracy when a computer solution 
was attempted. 

Following this degenerate result, the sampling geometry was altered by slight 
displacement of the locations of six of the frequency samples. The new sampling 
geometry is illustrated in Figure 3.4, with the circled points indicating the sample 
locations which differ from the degenerate case. The resulting filter frequency 
response is shown in Figure 3.5. 




Figure 3.4. Arrangement of Samples for Lowpass Filter 
Design, Case 2. 
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Figure 3.5. Frequency Response of Lowpass Filter 
Design, Case 2. 

Although the result may seem to be perfectly acceptable, the system of equa- 
tions involved in the solution can be considered to be ill-conditioned. Slight per- 
turbation of the position of a single sample point, as shown by the circled point in 
Figure 3.6, results in the unacceptable, excessively rippled response indicated by 
Figure 3.7. 

The above example demonstrates a negative aspect of the arbitrary-sampling 
approach. Although degeneracies, as defined, may seldom occur, ill-conditioned 
cases in which the determinant of matrix B in Eq. (3.12) is very small , but non- 
zero, may arise much more frequently. While specifying the frequencies u>i and u> 2 
at each sample point with a high degree of precision (analogous to placing samples 
on a very dense grid) may assure a very minimal probability of a degeneracy, it 
does not guarantee that the resulting system will not be ill-conditioned. 
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Figure 3.6. Arrangement of Samples for Lowpass Filter 
Design, Case 3. 



C. COMPUTATIONS 



A further negative aspect of this arbitrary-sampling method is the number 
of computations involved in determining filter parameters. While the uniform- 
sampling method of Chapter II was shown to be reasonably efficient in computation, 
the arbitrary-sampling method involves the solution of a rather large system of 
linear equations. This solution requires the inversion of an R x R matrix. This 
inversion dominates the computation process. It is well known that inverting an 
N xN matrix requires on the order of N 3 floating point operations, where a floating 
point operation consists of one multiplication and one addition. Since, for a filter 
with an N x x N 2 impulse response region of support, 



R - (M, + l)(A/ a + 1) - (— ^ 1- l)(— f l) 

= i(JV,+l)(JV 2 + !) 



(3.18) 
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Figure 3.7. Frequency Response of Lowpass Filter 
Design, Case 3. 



the order of the number of floating point operations involved in the arbitrary- 
sampling design process can be represented as 



C . ~ — N 3 /V 3 

^arb — 1 iV 2 * 

64 



(3.19) 



For the case in which = N 2 = N, this reduces to the following: 



c„„ = ^N‘. 

64 



(3.20) 



Thus, the arbitrary sampling method represents an approach which, while 
quite flexible, is computationally intensive and fails to provide guarantee of a 
unique solution. In the next chapter, an alternative nonuniform sampling approach 
is presented which, by placing some mild constraints on frequency sample location, 
always provides a unique design solution and is much more computationally effi- 
cient. 
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IV. A NOVEL 2-D NONUNIFORM FREQUENCY-SAMPLING APPROACH 



Chapters II and III presented two fundamentally different approaches to 2-D 
FIR frequency-sampling filter design. The uniform-sampling method possessed no 
flexibility in the location of frequency samples, yet was computationally efficient in 
determining filter parameters. The arbitrary sampling method, however, had few 
constraints on sample location but was computationally intensive. The approach 
described in this chapter represents a compromise. It is more constrained than the 
arbitrary sampling method, but has much more freedom in sample placement than 
for the uniform-sampling method. While not as efficient as the IDFT approach, 
this method requires significantly fewer computations when compared with the 
arbitrary sampling method. Additionally, existence and uniqueness of a solution 
for the impulse response is guaranteed, provided simple constraints on sample 
locations are met. 

A. APPROACH 

For convenience, the design method introduced in this chapter will be referred 
to as the “nonuniform sampling” method, although the arbitrary sampling method 
uses nonuniformly spaced samples as well. This method borrows ideas from the 
other two approaches. As in the arbitrary-sampling case, quadrant symmetry is 
specified. As was shown in Chapter III, the frequency response of a 2-D zero-phase 
FIR filter with an Ni X N 2 impulse response region of support can be represented 
as 

Af , M, 

H(ui,ui 2 ) = ^ a(n 1 ,n 2 )cosw 1 n 1 cosw 2 n 2 , (4-1) 

fli = 0 fl 3 = 0 
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where 



{ h(n i ,n 2 ) n k = 0, n 2 — 0 

2h(n 1 , n 2 ) / 0, ra 2 = 0 

2h(n 1 ,n 2 ) n x = 0, n 2 ^ 0 
4/i(n 1? n 2 ) ^ 0,n 2 0 



(4.2) 



with Mi = (iVi — l)/2 and M 2 = (N 2 — l)/2. In the arbitrary sampling method, 
the sequence a(n 1 ,n 2 ) was represented as a column vector, the basis functions 
rpi( w l5 w 2 ) were specified, and a system of equations of the form of Eq. (4.1) was 
solved directly. In the nonuniform sampling method, constraints placed on the 
sampling geometry result in Eq. (4.1) becoming a separable system. The system 
is then solved by a two-step process analogous to the row-column decomposition 
used in the uniform sampling method of Chapter II. The resulting systems of linear 
equations contain one-dimensional basis functions only; hence, linear dependencies 
associated with two-dimensional basis functions are not of concern. This results in 
assurance of a unique solution to the filter design problem. 

As in the arbitrary sampling method, samples are to be placed in region K = 
{(wi,w 2 ) : 0 < u>i < 7r;0 < ui 2 < ?r} of the two-dimensional frequency plane. 
As mentioned above, certain constraints on the location of frequency samples are 
required. Suppose any (M x + 1) distinct values of u> 1 are chosen on the interval 
0 < cui < 7 T. These values are denoted as {w lfc : k = 0, 1 , . . . , M x }. Then for each 
u lk , suppose any (M 2 + 1) distinct values of u 2 are chosen such that 0 < u> 2 < it. 
These values are denoted as {w 2< ( k ) : / = 0, 1, . . . , M 2 }. Then a unique solution for 
a(n i ,n 2 ) can be obtained from 

Mi M-, 

H{u lk ,u 2l (k)) = E E a(n 1 ,n 2 ) cosu^n! cos (k)n 2 ; 

n l = 0 n 3 = 0 



(4.3) 



0 <k<M 1 1 0 <1 <M 2 . 



Such a sampling geometry is shown in Figure 4.1 for the case M x = M 2 = 2. 
Notice that the selection of sample points is quite flexible. The values w 10 , Wn, 
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and cj 12 can be any arbitrary, distinct samples of uq in the region. Additionally, 
u> 20 (A:), w 21 (&), and w 22 (A:) represent any arbitrary', distinct samples of u> 2 for each 
k, k = 0,1,2. This is far less constrained than placing the samples on a uniform 
Cartesian grid. 




Figure 4.1. Nonuniform Sampling Example. 



Given the arrangement of samples described above, the sequence a(n 1 , n 2 ) can 
be found as follows. First, Eq. (4.1) is arranged as 

w, m , 



H( wi,w 2 ) = E[E a ( n l , n 2 ) COSWxUi J 

n 3 = 0 n i = 0 

A new function g(u; 1 ,n 2 ) is defined as 

M , 



COS U> 2 n 2 . 



(4.4) 



g(wi,n 2 ) = E a(n l , n 2 ) cos u) x n l . 



(4.5) 



ni = 0 
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The frequency response can then be expressed as 



Mi 



H((jJi,u 2 ) = Y: g(oJi,n 2 ) cos (jj 2 n 2 . 



(4.6) 



n ^ = 0 



For a particular value of namely ui lk , the response can be written as 



M 4 



ff(wik,w 2 ) = ^ g(u> lfc ,n 2 ) cosu> 2 n 2 . 



(4.7) 



n 2 = 0 



Now, for each w 2( (A;), l = 0, 1, . . . , M 2 , the above expression becomes 



M, 



II(uj lk ,uj 2l (k)) = ^ g{vik,n 2 )cosuj 2l (k)n 2 . 



(4.8) 



n q = 0 



Thus, the first step in the row-column decomposition is to solve a system of (M 2 +l) 
linear equations of the form of Eq. (4.8). In matrix form, this becomes 



H (wu,w 20 (k)) 
H{uj ifc,w 21 (k)) 




'1 cos u> 20 (k) 
1 cosu 21 (k) 


cos M 2 u> 20 (k) ' 
cosM 2 w 2 i(fc) 




ff ( v lk , 0 ) 

^(Wifc,l) 


- H(u lk ,u 2M , (k)) . 




.1 cosw 2M ,(fc) 


... cosM 2 w 2Wj (k ) . 




- ?( W lk > M 2 ) _ 



In vector-matrix notation, Eq. (4.9) is expressed as 

H fc = Ag fc , 



(4.9) 



(4.10) 



where H*. and g k are (M 2 + 1) x 1 column vectors, and A is an ( M 2 + 1) x (M 2 +1) 
square matrix. The vector g*. is then computed using the relation 



gfc = A 



(4.11) 



This entire process is performed for each k, k = 0,1,... ,Mi . Therefore, this step 
involves a total of (Mi + l) systems of (M 2 + l) linear equations each. 

Once the column vectors g k have been computed for all values of k, a matrix 
G of size (M 2 + l) x (M x + l) can be formed such that 

G = [ go gi ... i ] • (4-12) 
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Now the second step of the row-column decomposition can be performed. From 
Eq. (4.5), for a particular c u lk , 

M j 

g(u lk ,n 2 ) = ^ a(rii,n 2 ) cos (4.13) 

n i = 0 

For a particular n 2 , a system of linear equations can be formed by evaluating Eq. 
(4.13) for each k, k = 0, 1, . . . ,M X . In matrix form, this becomes 



' g(uio,n 2 ) 
g(oJn,n 2 ) 




"1 COSU7 10 

1 COS OJi 1 


cosMiCUio 

cosM^n 




a(0,n 2 ) ' 
a(l,n 2 ) 


-g(v 




.1 cosw 1Ml 


... cos Mi<jJi Ul . 




,a(Mi,n 2 ). 



Equation (4.14) can be expressed in vector-matrix notation as 

= Ba, (4.15) 

where is a (M x + 1) x 1 column vector formed from the elements of the n 2 th 
row of G, B is a (Mi + l) x (M x + 1) square matrix and a is a (M x + 1) x 1 column 
vector. Thus a is found from the relation 

a = B" 1 r n ,. (4.16) 

The vector a represents the values of a(n x , n 2 ) for a specified value of n 2 . To obtain 
the entire (Mi + 1) x (M 2 + 1) sequence a(nx,n 2 ), the process described by Eq. 
(4.16) is repeated for each value of n 2 , ra 2 = 0 , 1 ,..., M 2 . Therefore, a total of 
(M 2 + l) systems containing (M x + 1) linear equations each are involved in the 
second step of the method. Finally, the impulse response of the resulting filter can 
be evaluated using the relationship of Eq. (4.2). 

The following design example illustrates how this method can be employed. 
For this example, a 2-D circularly symmetric lowpass filter characteristic is desired 
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with unity gain in the passband region {(uq , uq) :_0 < \A-’i + u>\ < 0.27 r} and zero 
gain in the stopband region {(uq,u> 2 ) • 0.47T < y/wj + w 2 < 7r}. The region of 
support for the impulse response is specified to be 17 x 17. The arrangement of 
frequency samples shown in Figure 4.2 is to be used. Note that the sampling is 
not uniform. When samples are placed in the transition region, a linear slope 
characteristic between the two bands is used in specifying the desired frequency 
response value at those samples, as was done in the uniform sampling approach. 
Contour and perspective plots of the frequency response for the resulting filter are 
shown in Figure 4.3. 




Figure 4.2. Arrangement of Samples for 17 x 17 Lowpass Filter. 

One rather obvious extension of the nonuniform sampling method involves 
again specifying an Ni x N 2 impulse response region of support, but in this case 
performing the interpolation in the first step of the process on uq . Suppose any 
(M 2 + 1) distinct values of u> 2 are selected on the interval 0 < u> 2 < 7T. These are 
denoted as {w 2( , l = 0,1,..., M 2 }. For each value of w 2/ , any (Mi + l) distinct 
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(a) 




(b) 

Figure 4.3. Frequency Response of 17 x 17 Lowpass Filter. 
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values of w x are chosen on the interval 0 < uq <_7r. These values are denoted as 
{ uj i k (/), k = 0, 1, ... , Mi }. Such an arrangement of selected samples is illustrated 
in Figure 4.4 for the case M x = M 2 = 2. Now suppose /f(w n w 2 ) is specified at 
each selected sample {u xk (l),u 2 i) for k = 0, 1 , . . . , A/ x and / = 0, 1 , . . . , M 2 . 




Figure 4.4. Alternate Nonuniform Sampling Approach. 

Then there is a unique solution for a(n 1 ,n 2 ) which can be obtained from 

A i j A/ 3 

H(u lk (l),uj 2l ) = y2 a(n 1 ,n 2 )cosw lfc (/)n 1 cosw 2 ,n 2 ; 

m = o r»j = o ( 4 - 17 ) 

0 < k < Mi, 0 <1 <M 2 . 

The approach to the solution of the above problem and proof of its existence and 
uniqueness can be performed in a straightforward manner analogous to that of the 
basic approach. 
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It is worthy to note that the sampling geometry required by the uniform 
sampling method of Chapter II is, in reality, a special case of the general nonuni- 
form sampling form discussed here. Therefore, the results of the uniform sampling 
method can always be achieved by constraining the location of samples to the uni- 
form Cartesian grid. However, the method of solution is less efficient than the 
IDFT approach and so would probably not be used for this particular case. 

B. EXISTENCE AND UNIQUENESS 

A significant aspect of the nonuniform sampling method is that a unique solu- 
tion matching the specified frequency response at each sample point is guaranteed, 
provided samples are located as described. Proof of existence and uniqueness of a 
sequence h(n 1 ,n 2 ) which satisfies the interpolating conditions follows. 

The first step in the nonuniform sampling approach involves solution of (M y + 
1) systems of linear equations of the form of Eq. (4.9) or Eq. (4.10). A unique 
solution exists for each of these systems. Recall that Eq. (4.9) resulted from 
applying the interpolating conditions to Eq. (4.6). Through use of the Chebyshev 
polynomials, Eq. (4.6) can be expressed as 

M, 

H{u i,w 2 )= c ( w i, n 2 ) (cos w 2 ) n *, (4.18) 

n 2 = 0 

where a unique relation exists between the coefficients <7(ccq,n 2 ) and c{u> 1 ,n 2 ). 
Through the subtitution x = cos w 2 , the above becomes 

Ml 

if(w,,x) = c(u > 1 , n 2 )x”’ . (4-19) 

n ^ = 0 

The set {l,x,x 2 , . is known to be unisolvent on any interval [a,b] [Ref. 9: 

p. 31]. Therefore, a unique interpolation polynomial of the form of Eq. (4.19) can 
be found given any (M 2 + 1) distinct values of x. From the relation x, = cos w 2 , , 



45 



as long as distinct u> 2l (k) samples are specified on the interval {0 < ur 2 < 7r}, the 
corresponding x L will be distinct. Due to the unique relationship between Eqs. 
(4.6), (4.18), and (4.19), the interpolation implied by Eq. (4.10) will yield a unique 
solution for g fc , k = 0 , 1 ,..., M x . This implies that the matrix G defined in Eq. 
(4.12) will be unique. 

The second step of the nonuniform sampling approach involves solution of 
( M 2 + l) systems of linear equations of the form of Eq. (4.15). Each row of matrix 
B in Eq. (4.15) forms a vector of Chebyshev polynomials which form a unisolvent 
set. Therefore, B is invertible. Since a unique matrix G was found in the first step, 
the column vectors T n , are unique for each n 2 = 0, 1, . . . , M 2 . Therefore, a unique 
solution for a(n x ,n 2 ) exists which satisfies the overall interpolation problem. From 
the relationship indicated by Eq. (4.2), this implies existence and uniqueness for 
the resulting impulse response h(n 1 ,n 2 ). 

Ill-conditioned results normally occur in interpolation problems involving one 
variable when one sample is located in close proximity to another. Since the 
nonuniform sampling method uses a two-step process involving one-dimensional 
interpolating functions, the same result can be expected to hold here. Provided 
samples are not placed too close together, the resulting filter design should not be 
ill-conditioned. 

C. COMPUTATIONS 

While the nonuniform sampling method retains less flexibility than the ar- 
bitrary sampling method of Chapter III, it is more computationally efficient for 
filters of the same order. As in the analysis for the arbitrary sampling case, in- 
verting an N x N matrix is assumed to require on the order of N 3 floating point 
operations and additions, and inversion of the matrix is assumed to dominate the 
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